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Abstract 

We investigate the Coulomb dissociation of on ^osp]-, target at the beam 
energy of 250 MeV/nucleon, employing the cross sections for the radiative 
capture reaction ^Be(p,7)*B calculated within the Shell Model Embedded 
in the Continuum (SMEC) approach. In contrast to the situation at lower 
beam energies, the Coulomb breakup cross sections are found to be sensitive 
to the Ml transitions. Comparisons of SMEC and single-particle potential 
model predictions show that the Coulomb breakup cross sections at these 
high energies are sensitive to the structure model of '^B. Comparison with the 
preliminary data on the angle integrated spectra reported recently by the GSI 
group shows that the theory is able to reproduce the absolute magnitude as 
well as the shape of this data well. The contributions of the E2 component 
strongly depends on the range of the angles of the p — "^Be center of mass 
with respect to the beam direction, included in the angular integrations of 
the double differential cross section. If integrations are done over the angular 
range of 0°-1.87°, the E2 multipolarity contributes up to 25 % to the cross 
sections even for the relative energies of [p - ■^Be] below 0.25 MeV. However, 
these contributions are reduced by an order of magnitude at lower relative 
energies if the maximum of the angle integration is ~ 1°. 

PACS No. 25.70.De, 25.40.Lw, 96.60.Kx 



1 



Typeset using REVT^ 



I. INTRODUCTION 



The isotope produced in the Sun via the radiative capture reaction ''Be(p,7)®B is the 
principal source of the high energy neutrinos detected in the Super-Kamiokande and ^^Cl 
detectors In fact the calculated rate of events in the former detector (and also in the 
SNO experiment 0) is directly proportional to the rate of this reaction measured in the 
laboratory at low energies (~ 20 keV) 0. Unfortunately, the measured cross sections (at 
relative energies (Ecm) of [p — ^Be] > 200 keV) disagree in absolute magnitude and the 
value extracted by extrapolating the data in the region of 20 keV differ from each other by 
30-40 %. This makes the rate of the reaction '^Be(p, 7)^B most poorly known quantity in 
the entire nucleosynthesis chain leading to the formation of ^B . 

The Coulomb dissociation (CD) method provides an alternative indirect way to deter- 
mine the cross sections for the radiative capture reactions at low energies In this 
method, the radiative capture is reversed by the dissociation of the projectile (the fused 
system) in the Coulomb field of the target by assuming that the strong interaction between 
the nuclei is absent and the electromagnetic excitation process is dominated by a single 
multipolarity (see e.g. Refs. [§,|1)- However, in the CD of ^B, the contributions of E2 
and Ml multipolarities as well as nuclear breakup can be disproportionately enhanced in 
certain kinematical regimes pi]|JTT[] and a careful investigation is necessary to isolate 
the conditions in which these terms have negligible effect on the calculated breakup cross 
sections. 

For the CD measurements of ^B performed at RIKEN HTSHT^ at around 50 MeV/nucleon, 
it was found in a theoretical model ^ that the data are free from nuclear and E2 contri- 
butions if measurements are limited to Ecm < 0.70 MeV and to angles oi p—^Be cm. with 
respect to the beam direction {6cm) <4°. On the other hand, for measurements at sub- 
Coulomb energies [Q, these contributions turned out to be dominant everywhere P,p!7 . 



Furthermore, the multi-step breakup effects were also found to be quite important at these 
energies 

Recently, the Coulomb dissociation of ^B on a ^°^Pb target has been performed at GSI- 
Darmstadt at the beam energy of 250 MeV/nucleon |TP| , pn| . The advantages of the CD 



process at this high energy are: (i) measurements at Eqm lower than RIKEN are possible due 
to certain experimental advantages |T3J20]. (ii) the multi-step breakup processes {e.g. the 



Coulomb post-acceleration) are negligible ||T8|j21[| , and (iii) the nuclear breakup processes are 
expected to be negligible in comparison to the Coulomb one due to the strong enhancement 
in the virtual photon spectrum an advantage that existed already in the experiments 

performed at the RIKEN energies in the kinematical regime as mentioned above @ . At the 
beam energy of 250 MeV/nucleon, the E2 component is expected to be about one third of 
that found at RIKEN energies although, at the same time, the El component is also reduced 
by about one half. But more importantly, the Ml component is expected to be enhanced 
particularly in the first resonance region (0.5 MeV < Eqm < 0.65 MeV). This provides an 
alternative opportunity to test various models of the '^B structure, as they differ in their 
treatment of the resonance structure of this nucleus p2|-p5|]. 

The aim of this paper is to perform calculations of the CD of ^B on ^osp]^ target at the 
beam energy of 250 MeV/nucleon, using the cross section for the cature reaction '^Be(p,7)'^B 
calculated within a recently developed realistic shell model which includes coupling between 
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many particle quasi-bound states and the continuum of one particle scattering states (to 
be referred as the Shell Model Embedded into the Continuum (SMEC) approach) p6| , ^ . 
In order to investigate the sensitivity of the CD process to the structure of we also 
perform calculations within a single particle potential model using the parameters ginen 
by Esbensen and Bertsch (EE) I^Sj. We compare the results of our calculations with the 
preminary data on the angle integrated cross sections for this reaction taken recently in an 
experiment performed at the GSI, Darmstadt [ p]U| , pU[ . We discuss the kinematical regime 



where the breakup measurements performed at these energies can be reliably used to extract 
the astrophysical ^-factor (S'17) for the '^Be(p, 7)^B capture reaction. 

The remainder of this paper is organized in the following way. The formalism of the 
SMEC approach is described in the next section. The formulas used in the calculation of 
the Coulomb dissociation cross sections are also described here. The results of our numer- 
ical calculations and their discussions are presented in Sect. 3, while the summary and 
conclusions of our work is given in Sect. 4. 



II. THE FORMALISM 
A. The shell model embedded in the continuum (SMEC) 

In the next chapter, we shall study the CD process ®B — > ''Be-l-p, using the Shell Model 
Embedded in the Continuum (SMEC) which has been applied recently to examine structure 
for mirror nuclei, ^B, ^Li , and capture cross sections for mirror reactions ''Be(p, 7)^B, 
^Li(n, 7)^11 m,!?!. The SMEC model, in which reahstic Shell Model (SM) solutions for 
(quasi-)bound states are coupled to the one-particle scattering continuum, is a development 
of the Continuum Shell Model (CSM) []30H3^ for the description of low excitation energy 
properties of weakly bound nuclei. In the SMEC, the bound (interior) states together with 
its environment of asymptotic scattering channels form a quantum closed system. Using the 
projection operator technique, one separates the P subspace of asymptotic channels from the 
Q subspace of many-body states which are build up by the bound s.p. wave functions and by 
the s.p. resonance wave functions. P subspace is assumed to contain (A^ — 1) - particle states 
with nucleons on bound s.p. orbits and one nucleon in the scattering state. Also the s.p. 
resonance wave functions outside of the cutoff radius Rcut are included in the P subspace. 
The resonance wave functions for r < Rcut , are included in the Q subspace. The wave 
functions in Q and P are then properly renormalized in order to ensure the orthogonality 
of wave functions in both subspaces. 

For the (quasi-) bound many-body states in Q subspace one solves the SM problem 

Hqq^, = (1) 

where Hqq = QHQ is the SM effective Hamiltonian. What should be taken for coupling 
between bound and scattering states is in principle not known and we have decided to use 



a schematic combination of Wigner and Bartlett forces ||2^ 

Vu = -VS\a + (1 - a)Pr2]^(ri - ra) (2) 
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where is the spin exchange term and (1 — a) is the spin exchange parameter. We assume 
that effects of this couphng on the effective interaction in Q subspace are already contained 
in Hqq. 

The SM wave function has an incorrect asymptotic radial behavior for unbound states. 
Therefore, to generate both the radial s.p. wave functions in the Q subspace and the scatter- 
ing wave functions in P subspace we use the average potential of Saxon- Woods (SW) type 
with the spin-orbit part included 



[/(r) = Vo/(r) + V5o(41-s; 



1 dp{r) 
r dr 



Vr 



C 



where /(r) is the spherical symmetric SW form factor 



p{r) 



1 + exp 



(3) 



(4) 



The Coulomb potential Vc is calculated for the uniformly charged sphere with radius Rq. 
For the continuum part, one solves the coupled channel equations 



(+) 



E 







(5) 



where index c denotes different channels and Hpp = PHP . The superscript (+) means 
that boundary conditions for outgoing scattering waves are used. The channel states are 
defined by couphng one nucleon in the scattering continuum to a many-body state of (A^ — 1) 
- nucleus given by the SM. The channel - channel coupling potential in (RI) is 



H^. = {T + U)S,,^ + V 



(6) 



where T is the kinetic-energy operator and v^, is the channel-channel coupling generated 
by the residual interaction. U in (P) consists of an 'initial guess' U{r) given by and the 
diagonal part of the coupling potential v^^ which depends on both the s.p. orbit (pij and 
the considered many-body state J'^. Obviously, this correction cannot be neglected when 
generating s.p. wave function (pij for a given J'". These s.p. wave functions define Q subspace 
and thus modify the diagonal part of the residual force. This implies a self-consistent 
iterative procedure, because the change of s.p. wave function changes the correction coming 
from the residual force, and so on. This procedure yields the new self- consistent potential 



(7) 



and consistent with it the new renormalized matrix elements of the coupling force. The 
parameters of the initial potential U{r) are chosen in such a way that U^^^\r) reproduces 
energies of experimental s.p. states, whenever their identification is possible. 

The third system of equations in SMEC consists of inhomogeneous coupled channel 
equations 



(+) 



Hpp)u. 



(+) 



HpQ^ 



(8) 



with the source term Wi which is primarily given by the structure of - particle SM wave 
function $j . The explicit form of this source was given in These equations define 
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functions uul'^^ , which describe the decay of quasi-bound state $j in the continuum. The 
source Wi couples the wave function of - nucleon locahzed states with (A^ — 1) - nucleon 
locahzed states plus one nucleon in the continuum. Form factor of the source term is given 
by the self-consistently determined s.p. wave functions. 

The full solution is expressed by three functions $j , and Ui p7tp0 | 



1 



+ ^($, + uj,) _ < I Hqp I > (9) 
where 

Hqq = Hqq + HqpC'-p^HpQ (10) 

is the effective SM Hamiltonian which includes the coupling to the continuum, and G^"* is 
the Green function for the motion of s.p. in the P subspace. Matrix Hqq is non-Hermitian 
(the complex, symmetric matrix) for energies above the particle emission threshold and 
Hermitian (real) for lower energies. The eigenvalues, Ei — ^iTi , are complex for decaying 
states and depend on the energy E of particle in the continuum. The energy and width of 



resonance states is determined by the condition Ei{E) = E |30]. Inserting them in (^), one 
obtains 

+ y Qi — ^ ^ <^i\H\CE> (11) 

^ 'E-E, + {z/2)T,, ^11^^ ^ ^ 

for the continuum many-body wave function projected on channel c , where 

+ E dE Ce'^^^y^ < Ce' I H I > (12) 

is the wave function of discrete state modified by the coupling to the continuum states. It 
should be stressed that the SMEC formalism is fully symmetric in treating the continuum 
and bound state parts of the solution, \E'|; represents the continuum state modified by 
the discrete states, and VLi represents the discrete state modified by the coupling to the 
continuum. 



B. SMEC wave functions for 

The SMEC results depend mainly on (i) the effective nucleon - nucleon interaction in Q 
subspace, (ii) the residual coupling of Q and P subspaces, (iii) the self-consistent average 
s.p. potential which generates the radial form factor for s.p. bound wave functions and s.p. 
resonances. Cohen - Kurath (CK) interaction P3| is used for the SM effective interaction 



in Q subspace. The freedom of choosing asymptotic conditions in solving eqs. (|g) means 
that the zero on the excitation energy scale can be fixed arbitrarily. We choose the zero 
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on this scale by requiring that the lowest resonance {J'" = if) with respect to the proton 
emission threshold has its experimental position. An essential element of SMEC approach 
is the construction of Q - subspace. This is achieved by an iterative procedure which yields 
the self-consistent s.p. potential depending on the s.p. wave function (pij, the total spin J 
of the iV-nucleon system as well as on the one-body matrix elements of (A^ — 1) - nucleon 
daughter system. The parameters of initial SW potentials for different contributions of spin 
exchange component in the residual interaction are summarized in Table I. 

The unique potential U{r) is used for the calculation of self-consistent potentials for all 
many-body states in '^B, and for both lp3/2 and lpi/2 proton s.p. states. For neutrons, there 
is no correction from the residual interaction, and the average s.p. potential is chosen such 
that it yields 1^3/2 and lpi/2 neutron orbits at — 13.02 MeV and — 11.16 MeV respectively 



The quadrupole moment < Q > of is a useful test of the SMEC wave function. We 
have calculated < Q > following the approach of Carchidi et al. [Q. For the residual force 
(I) with = 650 MeV-fm^ and for the spin exchange parameter (1 — a) = 0.05 one 
finds 1^, < g >= 6.99 e fm^ in the good agreement with the experimental value ||35|| , 
< Q >= 6.83 ± 0.21 e fm . This theoretical value has been obtained assuming the effective 
charges, Cp = 1.35 , = 0.35, and the SM spectroscopic factors for the CK interaction. 

SMEC results depend sensitively on very small number of parameters. Some of them, like 
the parameterization of the residual interaction which couples states in Q and P subspaces. 



has been established previously |2gj23]. The others, related to the energy of s.p. states which 
determine the radial wave function of many-body states, are bound by the SM spectroscopic 
factors and experimental binding energy in studied nuclei. The spectrum of ^B depends 
strongly on couplings to the ground state (g.s.) of ^Be but changes very little if also the 
couplings to the excited state 3/2~ of ''Be are taken into account Width for J'^ = if 



state depends sensitively on the proportion of direct and spin exchange terms (see Table II) 
in the residual interaction. 

Varying the parameter of the spin exchange component for a fixed intensity of the cou- 
pling, we came to the conclusion that most satisfactory description of experimental data are 
achieved for small contribution of the spin exchange part, i.e., approaching the limit of pure 
Wigner force |3^ . This finding is consistent with the results of SM which strongly suggest 



an approximate validity of SU (4) symmetry in p-shell nuclei [ ^ , ^^ -|39|| 

SM energy of the first J'^ = 3f level is too low as compared to the experimental value 
(see Table II). The coupling to the continuum cannot correct for this deficiency. The width 
of 3f state differs by at least a factor 5 from the experimental data and here, again, the 
agreement between experiment and calculations improves when (1 — a) —^0 . There are 
several reasons for this discrepancy. Firstly, SM with CK interaction is not well describing 
energy of this state and, as pointed above, the width of the many-body state depends on its 
excitation energy with respect to the particle emission threshold (see Table II). Secondly, the 
wave function of experimental 3f state is certainly overlapping with the cluster configuration 
[^He — ^Li], which cannot be adequately described inp-space SM calculations. Experimental 
3^ state lies above the threshold for three-particle decay, *B — > [^He — p — ^He] . This decay 
channel largely contributes to the 3^^ width, but cannot be accounted for in the approxima- 
tion of one-particle scattering continuum. We have to keep in mind these limitations when 
analyzing the Coulomb dissociation cross-section for ^B . 
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C. Radiative capture 



The calculation of the capture cross-section in the SMEC goes as follows. The initial 
wave function for [p (g) ^Be] system is 



(13) 



and the final wave function for '^B in the g.s. ( J'^ = 2+) is 



•^f r . -|(J;) 

[Y'' X xT X X'* 



(14) 



It and s denote the spin of target nucleus and incoming proton, respectively, ^ifjj'l^ is the 



J t 

coefficient of fractional parentage and m^^^-^ is the s.p. wave in the many-particle state Jj . 
These SMEC wave functions, \&i(r) , ^/(r) , are then used to calculate the transition 
amplitudes T^"^ and T^^^ for El , E2 and Ml transitions, respectively |2^J2^. The radiative 
capture cross section is 



hcj \hcj 2s + l 2It + 1 
X ^ I x^^'^^^ 1^ 



X 



(15) 



a 



R2 



P 



1 



75 \kV \hcj \hcj 2s + 1 2L + 1 



X 



X 



E2 1 2 



(16) 



where /i stands for the reduced mass of the system. 

The astrophysical S-factor (Sn) is related to the capture cross section by 

Si7{Ecm) = ct^HEcm) Ecm e[2-''(^o«)], 

where rjlEcAi) is the Coulomb parameter and ttL represents the multipolarity. 



(17) 



1. SMEC results for '^Be{p,-ffB 

In Fig. 1, we show the astrophysical S-factors (^'17) (as defined by Eq. (2.17)) of El 
(solid lines), E2 (dashed-dotted line) and Ml (dashed line) multipolarities as a function 
of cm. energy for four versions of SMEC used in the calculations performed in this paper 



(see also Fig. 7 of Ref. |^). The parts (a)-(d) show the results obtained with versions 
I, II, III, IV of the SMEC model respectively. Version I corresponds to the spin exchange 
parameter (1 — a) in the residual coupling (|^) equal 0.27. This is the standard value resulting 
from a fit to the giant dipole resonance in ^^O ||30|j41[| . In version II, (1 — a) = 0.05, 
that corresponds to an almost pure Wigner force limit for this coupling. Strength of the 
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residual interaction in versions I and II of the SMEC model is, vj^^ = 650 MeV-fm^ and 
the parameters of the initial SW potentials are given in Table I. Versions III and IV are 
the same as versions I and II respectively except that the resonant Ml contribution for the 
J'^ = 3^ state has been omitted. This resonance has much smaller width than seen in the 
data due mainly to the missing three-body final state in the continuum, as explained in 
Sect. II. B. The Ml contribution and particularly its resonant part, are strongly dependent 
on the spin exchange parameter (1 — a) in the residual coupling (^. For (1 — a) = 0.05 
the resonant part of Ml transitions yields 5'-'^^ = 20.52 eV-b at the if resonance energy. 
This value is proportional to the square of spectroscopic amplitude of p-states, which for 
the CK interaction is —0.352 and 0.567 for pi/2 and P3/2 respectively. Ratio of E2 and El 
contributions at the position of if resonance is 8.15-10"^ or 7.72-10"^ depending on whether 
(1 — a) equals 0.27 or 0.05. The experimentally deduced value for this ratio, 6.7ll'l'^-10~'^ 
is consistent with both values of (1 — a). The E2 contribution contains both resonant and 
non-resonant contributions. Their ratio at l^" resonance is 0.187 (for {1 — a) = 0.05), and 
the contribution from different initial states, 0^, 1"'",2"'',3^ and 4+ is 1.17-10"^ /ib, 7.2-10"^ 
/ib, 4.64-10"^ /ib, 6.01-10^^ /ib and 5.64-10"^ /ib, respectively. Importance of this quantity 



has been suggested recently by Barker ^ 



El component provides the main contribution to the total capture cross-section. The 
low energy behavior of the astrophysical factor Sn{E) can be approximated by, Sn{E) = 
Su{0) exp{aE + PE"^). In the range of cm. energies up to 100 keV, it yields S'i7(0) = 19.594 
eV-b, a = -1.544 MeV-\ f3 = 6.468 MeV'^ for (1 - a) = 0.05. The above value of 5i7(0) 
is close to the values reported by Filippone et al. @] and Hammache et al. |0. At the 



position of 1^^ resonance, the calculated 5*17 - factor (517 = 40.67 eV-b) is somewhat smaller 
than that reported by Filippone et al. . 



D. Single- particle model for 

For a comparison, we shall also study the CD of ^B using a simple s.p. description of 
loosely bound proton in ^B, using the potential parameters given by Esbensen and Bertsch 
(EB) p8[. This will be referred to as the EB potential model in the following. The S-factors 
obtained with this model are shown by lines with solid circles in part (c) of Fig. 1. In 
this model, both the g.s. 2f as well as the resonances l^*" and 3f are assumed to have the 
structure [''Be(3/2^) (S> 1^3/2], the spectroscopic amplitudes for these states are assumed 



to be equal 1. This assumption is perhaps questionable for the li and Sf resonances |^ 
(the experimental spectroscopic factor for 1^ resonance for the mirror state in ^Li is 0.48 
). No intrinsic excitations of the ^Be are allowed and the s.p. potential well is adjusted 
for each state separately to reproduce either the one-proton separation energy (for the g.s. 
2f) or the excitation energies (for the l^" and 3f resonances). This extreme s.p. description 
of the '^B yields the width of resonances appreciably bigger than the experimental value, in 
contrast to the SMEC, which yields too narrow width for those states. For example, the 1^ 
width in the EB potential model is 70 keV as compared to the measured value 37 ± 5 keV 
and the 3^^ width is 1740 keV instead of 350±40 keV found experimentally [||]. In SMEC 
widths of l^" and Sf states are 26 keV and 35 keV, and 16.5 keV and 13 keV for for (1 — a) = 
0.05 and 0.27 respectively |]27[ . 
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In the EB potential model, the ratio S^'^/S^^ at the if peak energy equals 9.5 ■ 10~^ 
and situates at the upper most limit of the experimental values deduced by Davids et al. 
Ratio of resonant and nonresonant E2 contributions at this energy is 0.305. The 



astrophysical factor S'i7(0) is ~ 18.4 eV-b and is close to the value in the SMEC approach. 



E. Coulomb dissociation cross section 

The double differential cross-section for the Coulomb excitation of from its g.s. to 
the continuum, with a definite multipolarity of order vrA is given by 

2^ 77^,...^n, ^7 (^7). (18) 



In Eq. (17) ^Sg* defines the direction of the cm. of the p— '^Be system (to be referred as ^B*) 
with respect to the beam direction. a!^^{E.y) is the cross-section for the photo-disintegration 
process 7-l-^B — i>'^Be -|- p, with photon energy E^, and multipolarity vr = E (electric) or M 
(magnetic), and A = 1, 2... (order), which is related to that of the radiative capture process 
''Be + p ^ + •y through the theorem of detailed balance. E^ is given by Eqm = E^ + Q, 
with Q = 0.137. In most cases, only one or two multipolarities dominate the radiative 
capture as well as the Coulomb dissociation cross sections. nTr\{E^) in Eq. (0) represents 
the number of equivalent (virtual) photons provided by the Coulomb field of the target to 
the projectile, which is calculated by the methods discussed in Ref. [0,^^]. 



III. RESULTS AND DISCUSSION 

In Fig. 2, we show the results of the calculations for the energy differential cross sections 
for the reaction ^B -|- ^°^Pb ^B* -|- ^°'^Pb at the beam energy of 250 MeV/nucleon, using 
the capture cross sections obtained with versions I, II, III, IV of SMEC. 

Results shown in Fig. 2 have been obtained by integrating Eq. ([T8|) ) for 9^^, angles 
from 0.01° - 1.87°, which is the range of the angle integration in the preliminary GSI data 



as reported in [|19| , P0| . The cross sections for El, E2 and Ml multipolarities are shown by 
dashed, dotted and dashed-dotted lines respectively, while the solid line shows their sum. 
We note that all the four models reproduce the preliminary GSI data (which has been taken 
from the thesis of Federick Boue []Ty[) well. However, in models I and II, a peak appears at 
EcM ~ 1.2 MeV. This is because these models include also the contributions from the 3"'" 
state in the ^B continuum which in SMEC occurs at about 1.2 MeV excitation energy. This 
gives rise to a strong resonance in the Ml capture cross sections which is reflected in the 
corresponding Coulomb dissociation results. In models III and IV, contributions from this 
state is not included. In each case, the E2 multipolarity contributes to the extent of about 
25%. The maximum in the experimental cross sections is around Eqm = 0.5 MeV. The 
calculations cannot reproduce this without strong contributions from the Ml multipolarity, 
as will be shown later on. 

In Fig. 3, we present the comparison of the Coulomb dissociation cross sections (for the 
same reaction as in Fig. 2), obtained with capture cross sections of versions III (a) and 
IV (b) of the SMEC model and the EB potential model. The resuhs for El, E2 and Ml 
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multipolarities for SMEC models and the EB potential model are shown by dashed-dotted 
and dotted lines respectively. Their sum is shown by the solid and dashed lines in the 
two cases. We note that the SMEC cross sections are in somewhat better agreement with 
the preliminary GSI data. As compared to SMEC, the EB calculations under-predict the 
experimental data for Ecm larger than about 0.8 MeV, while they over-predict it around 
the 1+ resonance region. This can be traced back to the differences in predictions of the 
two models for the direct capture cross sections of various multipolarities (see Fig. 1). We 
note that while the El cross sections obtained with SMEC (version III) and EB model are 
similar, the two differ in case of SMEC version IV. On the other hand, SMEC E2 cross 
sections are always larger than that of the EB potential model. At the same time, the Ml 
cross section of the latter are much larger and wider in width as compared to that of former. 

Therefore, the Coulomb dissociation data at these high energies seem to show sensitivity 
to the capture cross sections calculated within different models of ^B structure. Although 
these are only preliminary data, yet SMEC models appear to be in a somewhat better 
agreement with it. Unfortunately, it is difficult to commment on the the difference seen in 
the widths of the Ml resonance in the two models. This is because the experiments have a 
finite resolution (~ 110 keV at Eqm = 0.63 MeV |2y]) that may give a large aparent width 
to the cross sections near the resonance region. Furthermore, the data are presented in the 
larger energy bins. It would, be worthwhile to improve upon both these aspects in the future 
studies. However, it must be added here that the width of this resonance is tested basically 
with the direct radiative capture (p, 7)data. 

We next investigate the role of the Ml multipolarity in the high energy data of GSI. 
We would like to recall that at lower beam energies {e.g., the RIKEN experiments [^,0]), 
the contribution of this multipolarity was almost negligible. In Fig. 4, we show the CD 
predictions (obtained with SMEC, version III) for El, E2 and Ml components of the an- 
gular distributions for the *B + ^"'^Pb ^B* + ^°^Pb reaction at the beam energies of 51.9 
MeV/nucleon. For completeness sake we also show the experimental data of Kikuchi et 
al. II^]. Note that since these CD calculations have been done within a pure semi-classical 
theory |47|, the agreement with the data, beyond 4° is not good as compared to that seen 
in Ref. 0. As has been discussed in Ref. 0, the point like projectile approximation of the 
semi-classical theory breaks down at angles beyond this. Inclusion of finite-size effects of the 
projectile reduces the cross sections at larger angles which leads to a better agreement with 
the experimental data 0. It may remarked here that although, the calculations reported in 
Ref. are done within a quantum mechanical theory, the point like projectile approxima- 
tion is still made there. Any how, the purpose of this figure is more to show the contribution 
of the Ml multipolarity to the CD cross sections. As can be seen, the contribution of this 
multipolarity is negligible even in the energy bin 500-750 keV. 

On the other hand, one can see from Fig. 5, where we show the CD calculations for only 
El (dashed lines) and E2 (dotted lines) multipolarities and their sum (solid line) for the 
same reaction as in Fig. 2, that at 250 MeV/nucleon, it is not possible to explain the data in 
the region of Eqm between 500-750 keV without the contribution of the Ml multipolarity. 
In this figure the results obtained with versions III (a) and IV (b) of SMEC are shown. The 
El + E2 cross sections in version IV are somewhat larger as the El component in this case 
is bigger. This sensitivity of the higher energy breakup data to the Ml multipolarity makes 
it possible to use this to supplement the information on the continuum structure of ^B which 
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was not feasible by similar studies at lower beam energies. 

In Fig. 6, we show the CD angular distributions {6 in this figure corresponds to ^8^*) 
calculated with capture cross sections of SMEC versions III and IV for the same reaction as in 
Fig. 2. These results have been obtained by integrating Eq. (^) over Ecm between 0.1 MeV 
to 3.0 MeV. The contributions of El, E2 and Ml multipolarities are shown by by dashed, 
dotted and dashed-dotted curves respectively. Their sum is shown by the solid lines. We can 
see that E2 and Ml contributions start becoming important already from 1°. Therefore, the 
requirement of the CD method that for a reliable extraction of the astrophysical ^-factor the 
data should be dominated by the excitation of a single multipolarity {El in present case), 
is more likely to be fulfilled in the measurements at GSI energies if the angle Oq^^, is kept 
below 1°. 

This point is further emphasized in Fig. 7, where we show the energy distributions of 
the CD cross sections obtained by integrating Eq. (17) for 9^^^ up to 1.87 ° (angular range 
of the data used in this paper) (part (a) ), and that obtained by performing the integration 
up to 1° only (part (b) ). 

The El, E2 and Ml components are shown by dashed, dotted and dashed-dotted curves 
respectively, while their sum is shown by the solid curves. The cross sections are shown only 
up to Ecm of 0.5 MeV, which is the region of interest for the determination of the S-factor. 
It can be seen from part (a) that if 6^^, goes up to 1.87°, the E2 components are substantial 
(up to about 25%) even at Ecm below 0.25 MeV. However, if 6*83* is confined to angles 
below 1°, the contributions of the E2 component are almost an order of magnitude down 
in comparison to those of El for Ecm below 0.30 MeV. Therefore, this provides a better 
possibility of a reliable extraction of 5*17. 



IV. SUMMARY AND OUTLOOK 

In this paper, we used the cross sections for the radiative capture reaction ^Be(p,7)®B, 
calculated within the shell model embedded into the continuum approach for the structure 
of ^B, to study the Coulomb dissociation of ^B on a ^osp]^ target at the beam energy of 250 
MeV/nucleon. Cross sections obtained with four versions of SMEC were used. Calculations 
were also performed with the capture cross sections obtained in a single particle model using 
the potential parameters given by Esbensen and Bertsch |^8[. Comparison of the calculations 



were made with the preliminary data for this reaction taken at GSI, Darmstadt recently. 

The CD cross sections at these high energy were found to be sensitive to the nuclear 
structure model of ®B. In contrast to the CD data taken at lower beam energies, the Ml 
multipolarity is quite important at higher beam energies. It may therefore, be possible 
to supplement the information on the continuum structure of *B from the CD studies at 
higher energies. As far as the preliminary data are concerned, the fits obtained with SMEC 
approach are somewhat better than those with the EB potential model. We noted that if the 
angles of the center-of-mass of the outgoing [p — ''Be] pair with respect to the beam direction 
were taken up to 1.87°, the E2 component is quite large even at [p — ^Be] CM energies 
below 0.25 MeV. To minimize the contribution of the E2 multipolarity, this angle should be 
confined to values below 1°. This conclusion appears to be by and large independent of the 
nuclear structure model of ®B. 
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TABLES 



TABLE L Parameters of the initial potentials U (r) ^ used in the calculations of self-consistent 
potentials C/('*'^)(r) for two parameters (1 — a) of the spin exchange term in the residual interaction 
(|2|). All these potentials have the same parameters of radius Rq = 2.4 fm, surface diffuseness 
a = 0.52 fm, and spin-orbit coupling Vso = — 4MeV. In all cases, the strength of the residual 
interaction (|) is ^^2°^ = eSOMeV-fm^ [2]. 



System 


ep3/2 [MeV] 


1 — a 


Vo [MeV] 


[p (8) ^Be] 


-0.137 


0.27 
0.05 


-40.045 
-37.660 



TABLE IL The dependence of B spectra on the relative strengths of direct and spin exchange 
parts of the residual interaction (^). Only ground state of ^Be was taken into account in all 
couplings. The proton separation energy is adjusted in order to reproduce the energy of the lowest 
resonance state 1^. The entries in this table are labelled by the value of the spin exchange parameter 
(1 — a) of the residual force. Strength of the residual interaction (|2|) is, Vi2^ = 650MeV-fm^. The 
cut-off radius is Rcut = 5fm except for the pi/2 s.p. wave function in 1^ many body states, which 
is in the continuum at about 300 keV above the threshold and for which larger cut-off was used 
Rcut = 10 fm. The numbers in parentheses are the widths of 3^ state if this state would be placed 
at the experimental energy. All units are in keV. 



State 


SM 


(1-a) 


= 0.73 


(1-a) 


= 0.95 


experiment 






energy 


energy 


width 


energy 


width 


energy 


width 


2+ 


-446 


-356 




-320 




-137.5±1.0 




1+ 


637 


637 


16.5 


637 


25.9 


637±6 


37±5 


3+ 


1246 


1294 


13.1 


1275 


34.9 


2183±30 


350±40 








(25.2) 




(67.4) 
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FIG. 1. The astrophysical S-factors for the radiative capture reaction '^Be(p,7)®B calculated 
with versions I, II, III and IV of the SMEC (shown in parts (a), (b), (c) and (d) respectively). Solid, 
dashed-dotted, and dashed curves represent the contributions of El, E2 and Ml multipolarities 
respectively. The results obtained with a single particle model with potential parameters taken 



from [28| are also shown together with model III. The El, E2 and Ml cross sections in this case 
are shown by the solid, dashed and dotted lines with solid circles. 
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FIG. 2. Comparison of the calculated Coulomb dissociation cross sections {da / (IEcm) with the 
experimental data for the breakup of on ^*^^Pb target at 250 MeV/nucleon, as a function of the 
p — Be CM energy. The results calculated with four versions of the SMEC are shown (models 
I, II, III and IV in parts (a), (b), (c) and (d) respectively). Dashed, dotted and dashed-dotted 
curves represent the contributions of El, E2 and Ml multipolarities respectively, while their sum 
is shown by the solid line. The experimental data are taken from the Ref. ||l^. These results 
have been obtained by integrating the double differential cross sections for angles in the range of 
0.0°-1.87°. 
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FIG. 3. Comparison of the CD calculations (for the same reaction and procedure as in Fig. 1) 
performed with the capture cross sections obtained with versions III (part (a)) and IV (part (b)) of 
the SMEC with those given by the EB single particle model The El, E2 and Ml components 
of the SMEC and EB models are shown by dotted and dashed-dotted lines respectively, while the 
sum of these components are shown by solid (SMEC) and dashed (EB) lines. 
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B + ^°^Pb 5 1 .9 MeV/nucleon 



10 



10 



10 



£cM = 500-750 keV 




£cM= 1250-1500 keV 




FIG. 4. El (dashed line), E2 (dotted line) and Ml (dashed dotted line) components of the 
Coulomb dissociation cross section eda/dO (calculated with version III of SMEC) as a function 
of the scattering angle 6 of ^B* for the dissociation of on 208p|j target at the beam energy of 
51.9 MeV/nucleon. The solid line shows their sum. Results for the relative energy bins of (a) 
500-750 keV, (b) 1250-1500 keV, (c) 2000-2250 keV are shown, e is the detector efficiency. The 
experimental data and the detector efficiencies are taken from 
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°B + "Tb 250 MeV/nucleon 

200 I , 




Er^ [MeV] 



FIG. 5. Comparison of the El + E2 (solid lines) CD cross sections calculated with the capture 
cross sections of versions III (part (a)) and IV (part (b)) of SMEC with the experimental data for 
the same reaction as in Fig. 1. The individual El and E2 components are shown by dashed and 
dotted lines respectively. 
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FIG. 6. Angular distributions for the CD of on ^ospj-, 250 MeV/nucleon calculated with 
versions III (part (a)) and IV (part (b)) of SMEC. These results have been obtained by integrating 
the double differential cross sections over the CM energies Ecm between 100 keV - 3.0 MeV. 



21 



+ ^^^Pb 250 MeV/nucleon 



10^ - 



> 



4P 10 



Eg" 1 



.5 



10 





. (b) 


li 1 , 1 .f' 1 , 1 , 


/" 

f 

_ 1, 

■ 1/ 

i ■ 1 , 1 ,.' 1 , 1 , 



0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 

FIG. 7. The Coulomb dissociation cross sections, da / dEcM ■, obtianed by integrating the double 
differential cross sections in the range of 0.01°-1.87° (part (a)) and 0.01°-1.0° (part (b)). The El, 
E2 and Ml components are shown by dashed, dotted and dashed-dotted lines respectively while 
their sum is depicted by the solid lines. 
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